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Abstract 

The dynamics of domain walls in nltrathin ferromagnetic strips with 
perpendicular magnetic anisotropy is studied from both numerical and an¬ 
alytical micromagnetics. The influence of the interfacial Dzyaloshinskii- 
Moriya interaction associated to a bi-layer strip arrangement has been 
considered, giving rise to the formation of Dzyaloshinskii domain walls. 
Such walls possess under equilibrium conditions an inner magnetization 
structure dehned by a certain orientation angle that make them to be 
considered as intermediate configurations between Bloch and Neel walls. 
Two different dynamics are considered, a field-driven and a current-driven 
dynamics, in particular, the one promoted by the spin torque due to the 
Spin-Hall effect. Results show an inherent asymmetry associated with 
the rotation of the domain wall magnetization orientation before reaching 
the stationary regime, characterized by a constant terminal speed. For a 
certain initial DW magnetization orientation at rest, the rotation deter¬ 
mines whether the reorientation of the DW magnetization prior to reach 
stationary motion is smooth or abrupt. This asymmetry affects the DW 
motion that can even reverse for a short period of time. Additionally, it is 
found that the terminal speed in the case of the current-driven dynamics 
may depend on either the initial DW magnetization orientation at rest or 
the sign of the longitudinally injected current. 


1 Introduction 

Domain walls (DWs) in ferromagnetic materials are boundaries separating re¬ 
gions (domains) that are magnetized differently. The efficient displacement of 
domain wall (DWs) along thin ferromagnetic strips is a prerequisite condition 
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for the application of the DW-based spintronic devices. mm DWs in soft ferro¬ 
magnetic strips, with Permalloy being the most commonly used material, have 
been extensively analyzed during the last decades, both from theoretical and 
experimental points of view (see |3] for an extended review). Nowadays most 
of the interest is focused on ultrathin strips sandwiched between a nonmag¬ 
netic heavy metal and an insulator, which depict high perpendicular magne- 
tocristalline anisotropy (PMA) and where the current-induced DW motion is 
anomalously efficient. |3] Due to the narrow DW widths and lower threshold 
current densities for motion of DWs, these asymmetric PMA systems are promis¬ 
ing platforms for solid-state magnetic devices based on electrically manipulate 
DWs. 

Apart from their potential for technological applications, and due to the 
rich physics involved, the analysis of the static properties of DWs in asymmet¬ 
ric stacks with strong PMA, and their field- and current-driven DW dynamics, 
are both also interesting from a pure fundamental point of view. The high effi¬ 
ciency of the current-induced DW dynamics was initially attributed to a Rashba 
effective held that stabilizes Bloch DWs against deformation, permitting high¬ 
speed motion through conventional non-adiabatic spin-transfer torque (STT). 
[4] However, a number of recent hndings suggest that STT contributes negli¬ 
gibly to DW dynamics in these ultrathin structures and interfacial phenomena 
are instead responsible. [SlISllTli Spin-orbit coupling (SOC) is at the basis of 
several effects in these systems. Firstly, SOC between the ferromagnetic strip 
and the insulating oxide overlayer induces the strong PMA in the ferromagnetic 
strip. Besides, SOC at the interface between the ferromagnetic layer and the 
heavy metal underlayer may result in different key phenomena explaining the 
experimental observations. [aiiiiTiiH] Indeed, the spin Hall effect (SHE) in the ad¬ 
jacent heavy metal has emerged as a possible alternative mechanism to the STT. 
The SHE produces a spin current from charge scattering in the heavy metal, 
and the resulting spin accumulation at the heavy-metal/ferromagnet interface 
generates a Slonczeswski-like torque O [TOl [H] sufficiently strong to drive the 
DW motion. [Sl|6l[7l[8] However, the SHE-induced torque alone cannot directly 
drive the magnetostatically preferred Bloch DWs in these materials. [T^] The 
Dzyaloshinskii-Moriya interaction (DMI), arising from the SOC and asymmet¬ 
ric interfaces, determines the magnetization texture of the DWs mill mils]. 
The DMI is a form of magnetic exchange interaction in which adjacent magnetic 
moments prefer to align orthogonal to each other with a certain handedness, in 
contrast with the ferromagnetic exchange interaction in which the magnetic 
moments prefer to align parallel. The DMI provides the missing ingredient to 
explain the current-induced DW motion stabilizing Neel DWs with a built-in 
longitudinal chirality, such that the SHE alone drives them uniformly and with 
high efficiently. piTginiiiHiiiiiEo! 

Although significant advances have been achieved during the last three years 
in the analysis of these systems, there exist still some interesting effects which 
have not studied so far. In particular, in asymmetric PMA heterostructures 
with moderate DMI interaction, the internal magnetization of the DWs do not 
depict neither pure Bloch (perfectly aligned along the transverse in-plane axis) 
nor pure Neel configurations (perfectly aligned along the longitudinal in-plane 
axis), but it adopts an intermediate state between them. m Here we report 
an analytical and micromagnetic study which indicates that the internal mag¬ 
netization moment is indeed degenerated, i.e., under moderate DMI, two DW 
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configurations are energetically possible. While the longitudinal magnetization 
component is governed by the chirality of DMI, the transverse component is 
independent of this chirality, so that it can take either a positive or a nega¬ 
tive value. This degeneration in the equilibrium state results in asymmetric 
held-driven and current-driven DW dynamics, which are both analytically and 
numerically analyzed here in detail. 

This work is structured as follows. In section the one-dimensional model 
is used to explore the equilibrium DW configurations at rest with emphasis on 
describing in details the introduced DW degeneration. The field-driven and the 
current-driven DW dynamics is evaluated in section]^ one-dimensional predic¬ 
tions are compared to full micromagnetic simulations. The main conclusions of 
our study are discussed in section]^ 


2 One-dimensional model 


The one dimensional model (IDM), as it was originally meant, [5T1|^ considers 
the existence of magnetic domains within the medium, separated by domain 
walls (DWs), transitional areas which establish an interface between two neigh¬ 
bor domains. These DWs can be schematically represented by a certain surface 
defined implicitly by an analytical function ip {x^y, z^t) = 0.[23) In fact, the 
expression ip (cc, y, z, t) = ipQ defines a collection of surfaces that propagates at 


every point in a direction given by the unit vector 


■0=bo 


at a speed that 


can be calculated as 


dtip 

IVV-I 


b='0o 


Since the orientation of the magnetization at every point is expressed by 
means of a couple of angles which are, in general, dependent on the space coor¬ 
dinates and time, there must exist a couple of extremal values for both angles 
defining each magnetic domain. In the simplest case, the orientation of the 
magnetization propagates coherently, and the interface is determined by a con¬ 
stant value of a certain component of the magnetization, which is given by any 
intermediate value of either one or both orientation angles. The instantaneous 
position of the interface can be then connected with a certain variable defined 
all over the interface trajectory. 

In the case of large strips, a IDM can be derived from the Gilbert equation 
of the magnetization along with the commonly named Walker trial’s functions, 
after the application of variational principles. [24] Walker trial’s functions in¬ 
troduce the concepts of DW position, DW magnetization orientation and DW 
width A, the latter obtained in terms of exchange and anisotropy free energies. 
In the absence of any external torque, one pair of coupled equations can be 
derived: 


_ i ^ 7 o da 

A “ 2fj.oMsAd^' 

<? , r 70 da 

O' — O =- 

A 2y,oMs dq ’ 


( 1 ) 

( 2 ) 


where 70 = 2.21 x 10 ®^ and a correspond respectively to the gyromagnetic 
ratio of the free-electron and the Gilbert damping constant, /tq is the vacuum 
permeability, and Ms is the saturation magnetization of the medium. The vari¬ 
able q represents the instantaneous position of the DW along the strip, and d) 
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determines in the case of strips with PMA the homogeneous orientation of the 
in-plane component of the magnetization with respect to the longitudinal direc¬ 
tion, as it will be further revised. Finally, tr represents the free energy density 
per unit surface associated to the presence of the DW, which is obtained by in¬ 
tegrating the volume free energy density along the longitudinal direction. This 
free energy density per unit surface originally included exchange, anisotropy 
and magnetostatic interactions. However, previous equations can be straight¬ 
forwardly used in order to add Zeeman interactions by considering them as part 
of such energy density. Some upgrades can also be made to include other exter¬ 
nal torques.pinning, m or thermal effects. pS] Besides, the suitability of 
the use of the abovementioned trial’s functions for strips with interfacial DMI 
has been adequately stated by other authors, m that introduced the concept 
of Dzyaloshinskii domain walls (DDWs). In this way, DDWs are intermediate 
cases when the orientation of the in-plane component of the magnetization with 
respect to the longitudinal direction differs from the extreme situations known 
as Neel DWs or Bloch DWs, as it will be further detailed. 

In this work, previous equations have been tailored to study the dynamic 
of deterministic DDWs in perfect strips with PMA driven by both external 
perpendicular field and SHE, being then written as follows for an up to down 
transition of the magnetization: 

$ (l -f a^) = 7 o {h^ + cos -b 

-b Q! 7 o sin$ (i/if cos $ — FId) , (3) 

^ (l -b a^) = a7o [h^ -b ^Hsh cos - 

— 7 o sin$ (7F^ cosd) — iF^)). (4) 

The different Fl-values stand for both the effective fields equivalent to the in¬ 
ner interactions within the ferromagnet, and the external stimuli. In this way, 
Hk = Ms {Nx — Ny), where Mg is the saturation magnetization and and 
Ny are respectively the so-called demagnetizing terms along the longitudinal 
and transverse axes, whose difference is approximately proportional to the film 
thickness t.[55] Additionally, Hjj is proportional to the Dzyaloshinskii-Moriya 
parameter D as Hd = — a • Finally, is the applied out-of-plane field 
and Hsh defines the spin-orbit torque (SOT) associated to SHE, depending on 
the Spin-Hall angle in the heavy metal and being proportional to the 

longitudinally injected current as Hsh = with e representing the 

elementary charge, and h the reduced Planck constant. 

2.1 Equilibrium condition 

Previously to the dynamic characterization, a brief review of the equilibrium 
condition for DDWs is needed. In the absence of external stimuli the DW 
magnetization is oriented along a certain angle $1^^ = tpo- The equilibrium 
condition is then derived from the minimization of the free energy density per 
unit surface as calculated in m with regard to the orientation angle. This leads 
to the equations: 

sin (fo {Hk cos ifo - Hd) =0, 

Hk (cos^ (po — sin^ po) — Hd cos < 0, 
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(5) 

( 6 ) 




that admit a solution in the form: 


cos (po = sign 


Er) 

Hk) 


min 




(7) 


which includes the extreme cases of absence of DMI, that is, Ho = 0 and 


ipo = , as for Bloch DWs, or strong DMI, so that 


Hr 


Hk 


> 1, and ipo = 0 


or (^0 = as for Neel DWs. Any other values of the DMI lead to DDWs, so 
that two orientations of the magnetization within the wall are possible for every 
single value of the DMI parameter, being these two orientations symmetric with 
respect to the longitudinal axis. This is schematically depicted in figurej^ where 
some axes have been proposed in order to clarify the geometry and the mag¬ 
nitudes here defined. A strip with three micromagnetically calculated domains 
(up-down-up) is presented, with corresponding to the out-of-plane compo¬ 
nent of the magnetization, and rrix and my being respectively the longitudinal 
and transverse in-plane components. Details on how micromagnetic calcula¬ 
tions have been carried out are further given along the paper. The stacked 
color maps in each subfigure show the strength of any of these components at 
the domains and at the DWs. Since the out-of-plane component of the mag¬ 
netization goes through zero within the DW, the in-plane component take at 
such point a maximum of magnitude mow oriented along the $ angle. As it 
has been stated, under equilibrium conditions this angle is fixed to one of the 
two possible values of ipo which are a solution of eq. [Tj In this way, while the 
longitudinal component of the magnetization within the DW is governed by the 
chirality due to DMI, the transverse component can take two symmetric values. 
In the case of figure[^a) both DWs have a positive transverse component, while 
in figure [^c) both DWs have a negative transverse component. Different signs 
of the transverse component for neighbor DWs are also possible, as it is shown 
in figure [^b) 

Next section shows that these two degenerated equilibrium states of DDWs 
for a given value of the DMI parameter present dynamic behaviors under the 
same external stimulus, such as out-of plane applied fields or SHE associated 
to longitudinally injected currents, which are intrinsically different. These dif¬ 
ferences can be also revealed by applying the external stimulus in opposite 
directions to a certain DDW, as it is further demonstrated. 


3 Dynamics of Dzyaloshinskii DWs 

As the above referenced experimental results evidence, the dynamic of a DW 
admits a stationary motion under the application of certain stimuli, which is 
characterized by a constant DW speed. The range of the applied stimuli that 
allows such a stationary behavior depends on the type of torque or interaction. 
For example, field-driven dynamics is characterized by two different regimes, 
the abovementioned stationary one for low applied fields, and a lower-velocity 
regime with a region of negative mobility just above the critical field called 
Walker field. Nevertheless, current-driven dynamics as this induced by the SHE 
is characterized by a single and stationary regime. This can be straightforwardly 
deduced from eq. ([^, which states that <!> is a continuous and bounded function 
of d>. If the DW starts from an equilibrium condition, the immediate application 
of either external stimulus gives to $ a certain finite value, whose sign defines if 
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Figure 1: Definition of the geometry and some of the magnitudes involved in 
this study, in particular, the applied stimuli in the form of out-of-plane magnetic 
fields Hz and longitudinal currents jx- Some DDW static configurations are 
also depicted. In particular, subhgures (a), (b) and (c) show how the chiral 
character of the DMI forces the longitudinal component of the magnetization 
rux within a DW to rotate either clock- or counterclockwise, but not both. 
However, the transversal component of the magnetization rriy is achiral, so that, 
this component is free to rotate within the DW. In this way, this rotation can 
take either the same direction within two consecutive DWs, as in subfigure (b), 
or different directions, as in subfigures (a) and (c). 


the DW magnetization orientation initially rotates clock- or counterclockwise. 
Two behaviors may occur from this point on. On the one hand, $ may vary 
without a change of sign as $ varies in time. This situation occurs when both the 
upper and lower bounds of $ have the same sign, then limiting the rotation speed 
of the DW magnetization orientation, but preventing the stationary regime from 
being reached. That is the case of large applied out-of-plane helds, that may 
mask the other terms in the RHS of eq. ([^, and Walker breakdown occurs. On 
the other hand, the rotation of the DW magnetization orientation may reach a 
certain angle so that $ goes through zero. From this instant on, rotation stops 
and the stationary motion is reached. 

According to the previous discussion, it can be stated that the stationary 
motion is characterized through the conditions $ = 0, and then q = const, the 
former leading to the relationship: 

TT 

Hz + cos = a sin (ps (Hd - Hk cos ps ), (8) 

Ps representing the DW magnetization orientation when the stationary regime is 
reached. From the comparison between eqs. 0 and ([^ , and the characteristics 
and dependences of $ given by eq.(j^, it can be immediately inferred that the 
transition from the DW equilibrium state to its stationary motion involves a 
monotonous rotation of the DW magnetization orientation from p^ to ps- 

3.1 Field-driven dynamics 

As abovementioned, the DDW dynamics driven by the application of an out- 
of-plane external field admits two different regimes. For the sake of a lighter 
notation, let us rewrite eq.Q in the absence of SHE as: 

h = s\nps{5 — cosps), (9) 
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where h = is a normalized value of the applied field, and 5 = defines 
the normalized strength of DMI. Since the RHS of previous equation is bounded, 
there are h-values so that the equation has no solution. These /i-values are then 
above Walker breakdown, and the frontier can be delimited by considering the 
global maxima of eq. (§. At the frontier, the DW magnetization orientation tpvv 

= 0. This is a necessary, but not sufficient condition. 

Vb—<Pw 

In any case, this condition leads to the relationship: 


is such so that 


cos (fw = 




( 10 ) 


For strong DMI, characterized by absolute values of <5 greater than one, only 
one sign in the equation above leads to a valid result, i.e., the minus sign for 
positive (5, then defining the positive and negative Walker fields. It must be 
noted that the equilibrium condition for such i5-values establishes pure Neel 
walls, that present a perfectly symmetric behavior under the proposed stimuli, 
as it is further discussed. However, in the case of weak DMI, when DDWs are 


present, both signs in eq.(lO) are possible, one of them corresponds to global 


extrema while the other correspond to local extrema. As in the case of large 
<5, global extrema define the limit of the stationary motion under an external 
applied field, i.e.. Walker breakdown. 

In order to analyze the importance of local extrema, let us consider a DDW 
initially in equilibrium with a DW magnetization orientation fulfilling the con¬ 
dition cos:/5o = according to eq.Q and the definition of 5. The application 
of the external field h gives rise to a counterclockwise rotation for positive h, 
that is, applied field along the direction of the magnetization in the up-domain, 
and a clockwise rotation for negative /i, as it can be inferred from eq.(|^. For a 
certain sign of h, and provided h is sufficiently strong, the DW magnetization 
orientation can reach this local extremum. However, if this /i-value is fixed, but 
its sign is changed, no local extrema are reached. This is clearly shown in Figure 

Plots represent the DDW dynamics under the influence of an out-of-plane 
external field. In order to depict realistic values of time, DW position q and DW 
speed V, a value of A of about 6nm has been considered in the graphs. An initial 
DW magnetization orientation ipo = 60°, corresponding to a (5-value of 0.5, has 
been taken. Figure [^a) is obtained for applied fields of h = ±0.23, slightly 
higher in absolute terms than the value need to reach a local minimum. A clear 
asymmetry with the sign of h is shown. While the positive h-value promotes 
a soft transition from the equilibrium condition to the stationary motion, the 
negative h-value gives rise to an abrupt transition of the DW magnetization ori¬ 
entation prior to reach the stationary motion. During this transition, the DDW 
motion even reverses for a short period of time, as shown in the evolution of the 
DW position q. Similar results can be obtained if the applied field is increased 
close to the Walker breakdown (Figure [^b). However, the shown asymmetry 
do not affect to the terminal speed which is given by joH^A. 

Micromagnetic simulations have been carried out with the help of the GP- 
Magnet software packagein order to support these analytical results. Along 
these simulations, a ferromagnetic strip of Inm thickness and I60nm width with 
PMA grown on a heavy metal layer has been considered. Values of the material 
parameters that are usually found in the literature have been taken: saturation 
magnetization of Ms = 7 ■ 10® — , exchange constant of A = 10“^® — , uniaxial 
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Figure 2; DDW dynamics under the influence of an out-of-plane external field. 
In these graphs, h represents a normalized value as it has been defined along 
the text. Positive h-values stand for fields applied along the direction of the 
magnetization in the up-domain, while negative values stand for fields applied 
in the opposite direction, i.e., the magnetization in the down-domain, q rep¬ 
resents the instantaneous DDW position, and v is the instantaneous speed. <i> 
corresponds to the DW magnetization orientation. Plots correspond to the nu¬ 
merical calculation of eqs.Q and Figure a) is obtained for applied fields of 
h = ±0.23, and figure b) is obtained for applied fields of h = ±0.88. The sign 
of the applied field may promote completely different DDW dynamics. 































Figure 3: Comparison between micromagnetic simulations and the predictions of 
the IDM analytical model. Different signs of the applied fields lead to different 
behaviors of the a DDW. 

anisotropy constant ku = 4.8 • 10® A DMI parameter D of the order of 
10“®^ has been used. 

As an example, Figure [^presents the micromagnetic results obtained in the 
case of a sample with 5 = —0.5. In this case, equilibrium conditions determine 
that one of the two possible orientations at rest of the DW magnetization is 
tpo = 120°, as it can be checked in the graph for h = 0. The figure shows 
the dependence of the DW magnetization orientation at stationary motion as 
a function of the external field h. Dots are calculated by means of micromag¬ 
netic simulations, while the continuous plot corresponds to the IDM analytical 
treatment. In this example, the stationary motion is smoothly reached from 
the equilibrium orientation for negative helds, up to the corresponding negative 
Walker field. However, an abrupt transition occurs for positive h-values, before 
the Walker held is reached. The inset shows the normalized in-plane compo¬ 
nents of the magnetization at the center of the DDW. is the component 
along the longitudinal axis, while rriy is the transverse component. Starting 
from equilibrium conditions (hlled dot), negative helds give rise to a smooth 
clockwise rotation of the DW magnetization orientation for all applied helds up 
to the Walker held. On the contrary, positive helds leads to a counterclockwise 
rotation that pass through an abrupt reorientation of the DW magnetization 
prior to reach Walker breakdown. 

In order to complete the IDM analytical study, Figure shows the depen¬ 
dence of the DW magnetization orientation at stationary motion on the normal¬ 
ized out-of-plane held h with 6, i.e., as a parameter. On the one hand, large 
ratios promote Neel walls, their chirality depending on the sign of the DMI 
parameter D, whose behavior is completely symmetric with regard to the sign 
of h. As it is well-known, the larger the ratio is, the farther the two endpoints of 
the corresponding curve are from each other, meaning that Walker breakdown 
requires larger applied helds. On the other hand, no DMI allows the existence of 
Bloch walls. In this case, the dependence of the DW magnetization orientation 
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Figure 4: Dependence of the DW magnetization orientation at stationary motion 
on the applied out-of-plane field with as a parameter. While both Bloch 
and Neel walls present symmetric behavior, DDWs present a clear asymmetry 
with the applied field governed by the strength of the DMI. 


at stationary motion is depicted for a Bloch DW with tpo = 90°. The behavior 
is also symmetric with regard to the sign of h. However, DDWs present a clear 
asymmetry under the application of an external out-of plane field. This behav¬ 
ior can be considered as hysteretic in the following sense. Under equilibrium 
conditions, a DDW may take two possible orientations, symmetric with respect 
to the longitudinal axis. The application of an out-of-plane field promotes the 
rotation of the initial DW magnetization orientation before a state of station¬ 
ary motion is reached. Depending on the sign and the strength of the applied 
field, the DW magnetization orientation may undergo an abrupt reorientation 
prior to the stationary state. If this reorientation occurs, the DDW do not 
recover its original orientation at rest, but the symmetric one with respect to 
the longitudinal axis, after field removal. The existence or absence of these DW 
magnetization reorientation noticeably defines DDW dynamics, as the dynamic 
plots in Figure [^showed. 

3.2 Current-driven dynamics 

As it has been already stated, field-driven dynamics and current-driven dynam¬ 
ics governed by SHE present a noticeable difference, that is, the absence for 
the latter of Walker breakdown. This can be straightforwardly proved from 
eq.Q. Analogous to the precedent study, a lighter version of this equation can 
be proposed in the following way: 

/i = tan (/3s ((5 — cos :/3s), (11) 

where now h = ■ Since the RHS of previous equation is not bounded, 

a stationary motion is always possible under the SHE torque. However, the 
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Figure 5: Current-driven DDW dynamics under the influence of SHE for positive 
and negative longitudinal currents as calculated analytically with the help of the 
IDM. As in the case of the field-driven dynamics, asymmetry with the sign of 
the stimulus is also noticeable. 

function defined in eq. 0 may possess local extrema for absolute (5-values 
ranging from 0 to 1, that is, for DDWs. This local extrema are reached when 
the DW magnetization orientation fulfills the condition cos ipL = As in the 
precedent case, when the DW magnetization orientation reaches either of the 
two local extrema, an abrupt reorientation occurs prior to the establishment of 
the stationary motion. 

As a first example. Figure shows the DDW dynamics under the SHE for 
a certain current flowing longitudinally in both one and the opposite direction 
calculated by means of the IDM analytical description. In order to give a much 
more clear idea about orders of magnitude, one set of concrete values has been 
chosen, which are current density jx = ±4-10^° Spin-Hall angle 6sh = 0.11, 
with a damping constant a = 0.013 and a DW width A of about 6nm, values 
that can be also found in the referenced literature. In a similar fashion to the 
field-driven dynamics, a smooth rotation of the DW magnetization is promoted 
if the current flows in one direction, while an abrupt transition occurs when the 
current flows in the opposite direction. 

Differently from the field-driven dynamics, the IDM also predicts the asym¬ 
metry of the terminal speed at stationary motion, which undergoes a sharp tran¬ 
sition for currents that lead the DW magnetization orientation to values close to 
the local extrema. Figur^ shows the dependence of this terminal speed on the 
/i-value defined after eq.([H). Dots correspond to numerical simulations carried 
out with GPMagnet using the abovementioned parameters, while the continuous 
plot has been analytically calculated from the IDM. A rather good agreement 
between numerical and analytical calculations is obtained. Both show the sharp 
transition predicted by the IDM, in this case, for an /i-value around -20. It can 
be noticed that the terminal speed at the peak of this sharp variation, calcu¬ 
lated from the IDM as ^Hsh, may surpass the terminal speed for infinite 

current given by ^oHjjA. 

As a summary, the analytically calculated dependence of the DW magne¬ 
tization orientation at stationary motion on the applied current with as a 
parameter is depicted in Figure The case of Bloch DWs is not considered, 
since SHE do not promote their motion in the system under study. For 
ratios greater or equal to one, the plots show that Neel walls behave symmet¬ 
rically under the influence of the torque due to SHE. However, DDWs have a 
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Figure 6: Dependence of the DDW terminal speed at stationary motion on the 
applied longitudinal current due to SHE. The h-value in the graph is defined as 
proportional to the current density (see text). The graph compares micromag- 
netic simulations and the results obtained from the IDM model. A noticeable 
asymmetry is found if the stimulus is reversed, leading to a sharp transition of 
the terminal speed in a certain range of applied currents. 


hysteretic-like behavior, rather analogous to the one found for the field-driven 
dynamics. It must be noticed that, according to this hysteretic-like behavior, 
the terminal speed depends on the initial DW magnetization orientation under 
equilibrium conditions, i.e., for the same applied longitudinal current, the dy¬ 
namics depends on whether the initial DW magnetization orientation at rest is 
pointing in one of the two possible directions, symmetric with respect to the 
longitudinal axis, that this orientation can take. 

4 Conclusions 

The dynamics of Dzyaloshinskii DWs in ferromagnetic strips with PMA has 
been analytically studied with the help of the IDM, and the results compared 
with micromagnetic simulations. Two different stimuli have been considered: 
out-of-plane applied magnetic fields and torques induced by longitudinal cur¬ 
rents due to SHE. The application of such stimuli may lead to a stationary 
regime characterized by a constant terminal DW speed. However, the behavior 
of the DDW dynamics has been found to be asymmetric with respect to either 
the sign of the applied stimulus, or the initial DW magnetization orientation. 
This asymmetry arises from the fact that, depending on this sign, the DW mag¬ 
netization orientation may rotate clock- or counterclockwise. The stationary 
regime is then reached when this rotation stops. However, the transition from 
the DW magnetization orientation at rest to the orientation at stationary mo¬ 
tion can be either smooth or abrupt depending on the rotation. The initial DW 
magnetization orientation under equilibrium conditions then determines which 
sign of the stimulus promotes its abrupt reorientation prior to reach the sta¬ 
tionary regime. This process can be characterized in the presented results by a 
reverse of the DDW motion for a short period of time. 
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Figure 7: Dependence of the DW magnetization orientation at stationary motion 
on current through the normalized /i-value (see text) with as a parameter. 
It is shown that Neel walls present symmetric behavior, while DDWs present a 
hysteretic-like behavior. As a consequence, a DDW may reach different terminal 
speeds depending on either the sign of the current or its initial orientation at 
rest. 


Additionally, some differences between field-driven and current-driven dy¬ 
namics have been reported. Some of them have been already highlighted by 
other authors, such as the absence of Walker breakdown in the case of current- 
driven dynamics governed by SHE. Besides, the asymmetry here reported do 
not affect to the DW terminal speed in the case of field-driven dynamics. How¬ 
ever, the terminal speed for current-driven dynamics is found to be asymmetric 
in the range of lower injected currents, then depending on either the initial DW 
magnetization orientation or the sign of the applied stimulus. 

Finally, since for higher injected currents the abrupt reorientation of the 
DW magnetization occurs, the degenerated condition of the initial DW mag¬ 
netization orientation also determines different traveled distances for a train of 
current-driven DWs in a strip, which may be also of experimental relevance. 
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